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Abstract 

o 

We survey combinatorial interpretations of some dozen identities for tfie double 
I factorial such as, for instance, (2n — 2)!! + Ym=2 ^'^^~{2k-\)\\^^^'^' ~ (^'^ ~ ■*■)'" 

^ ' methods are mostly bijective. 



1^ ■ 1 Introduction 
m 

There are a surprisingly large number of identities for the odd double factorial (2n — 1)!! = 
^ ■ {2n — 1) ■ {2n — 3) ■ ■ - 3 • 1 = that involve round numbers (small prime factors), as 

well as several that don't. The purpose of this paper is to present (and in some cases, 
review) combinatorial interpretations of these identities. Section 2 reviews combinatorial 
' constructs counted by (2n — 1)!!. Section 3 uses Hafnians to establish one of these manifes- 
tations. The subsequent sections contain the main results and treat individual identities, 
presenting one or more combinatorial interpretations for each; Section I is devoted to 
round-number identities. Section to non-round identities. Section 6 to refinements in- 
volving double summations interpreted by two statistics in addition to size. 

The even double factorial is (2n)!! = 2n ■ {2n — 2) • ■ - 2 = 2"n!. The recurrence k\\ = 
k{k — 2)!! allows the definition of the double factorial to be extended to odd negative 
arguments. In particular, the values (—1)!! = 1 and (—3)!! = — 1 will arise in some of 
the identities. We use the notations [n] for the set {1,2, . . . ,n} and [a,b] for the closed 
interval of integers from a to b. For negative k, (2) = as usual, except that identities 
(2) and (10) below require (:J) := 1. It has become customary to draw trees down 
but when a construction involves growing a tree, it seems more natural to draw it up. 
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Also, to visualize ordered trees as Dyck paths, they must go up. So we will combine 
arborological pictures conjuring roots and leaves with the usual genealogical terminology 
of parents, children, and siblings. We sometimes refer to (clockwise) walkaround order of 
the edges/vertices in a tree; more formally, it is the order edges are visited in depth-first 
search and the preorder of the vertices. 

2 Combinatorial manifestations of (2n — 1)!! 

We begin with a review of some combinatorial manifestations of (2r2 — 1)!!, each illustrated 
for the case n = 2. In all cases, the parameter n is the size of the object. 

2.1 Trapezoidal words 

The Cartesian product [1] x [3] x . . . x [2n — 1]. 

11, 12, 13. 

The elements of this Cartesian product, the most obvious manifestation of {2n — 1)!!, 
were called trapezoidal words by Riordan [ ]. A minor variation is symmetric trapezoidal 
words: the Cartesian product [0, 0] x [—1, 1] x [—2, 2] x . . . x [— (n — 1), n — 1)]. 

2.2 Perfect matchings 

Perfect matchings of [2n]. 

12/34, 13/24, 14/23. 

A perfect matching of \2n\ = {1,2, ...,2n} is a partition of [2n] into 2-element sub- 
sets a(l) 6(1) / a{2) 6(2) / . . . / a{n) b{n) written, in standard form, so that a{i) < b{i) 
for all i, and a(l) < a(2) < . . . < a{n). Erasing the virgules (slashes) then gives a 
bijection to the perfect matching permutations of [2n], denoted P„: the permutations 
a(l) 6(1) a(2) 6(2) . . . a{n) h{n) (denoted (a, b)) of pri] satisfying a{i) < h{i) for all i, and 
a(l) < a{2) < ... < a{n). Given a perfect matching permutation (a, b), form a list 
whose a(i)-th and 6(i)-th entries are both i. This is a bijection to the permutations of 
the multiset {1, 1, 2, 2, . . . , n, n} in which the first occurrences of 1, 2, . . . , n occur in that 
order. 

A perfect matching of [2n] can be regarded as a fixed-point-free involution on [2n] 
and also as a 1-regular graph on [2n], whose pictorial representation is sometimes called 
a Brauer diagram. 
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2.3 Stirling permutations 

Permutations of the multiset {1, 1,2,2, ... ,n, n} in which, for each i, all entries between 
the two occurrences of i exceed i [2]. 

1122, 1221, 2211. 



2.4 Increasing ordered trees 

Increasing vertex-labeled ordered trees of n edges, label set [0,n] [3, 4]. 

2 

12 2 1 I 

V V I' 



Given such a tree, delete the root label and transfer the remaining labels from vertices 
to parent edges. Walk clockwise around the tree thereby traversing each edge twice and 
record labels in the order encountered. This is a bijection to Stirling permutations due 
to Svante Janson [ ]. Label sets other than [0,n] may arise and the term standard then 
emphasizes that a tree's label set is an initial segment of the nonnegative integers. 



2.5 Leaf-labeled 0-2 trees 

0-2 trees (rooted, unordered, each vertex has or 2 children) with n + 1 labeled leaves, 
label set [l,n+ 1] [ , Chapter 5.2.6]. 



12 13 2 3 




2.6 Height-labeled Dyck paths 

Dyck paths [ , Exercise 6.19] of n upsteps and n downsteps in which each upstep is labeled 
with a positive integer < the height of its top vertex. 




This item is due to Jean Frangon and Gerard Viennot [7, 8]; they observe that it is 
a consequence of their bijection from permutations to certain marked-up lattice paths. 
Here is perhaps the simplest proof and several further bijections that prove the result will 
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appear in the course of the paper. A height-labeled (HL) Dyck path P of size n — 1 has 
2n — l vertices, each of which can be used to construct a height-labeled Dyck path of size 
n: split P at the specified vertex into subpaths Pi and P2, insert an upstep between Pi 
and P2, increment by 1 the labels on P2, and append a downstep, as illustrated below. 




HL Dyck 5-path with marked vertex HL Dyck 6-path 



The process can be reversed by locating the last upstep with label 1 in a size-n path. 
Thus there is a multiplying factor of 2n — 1 from size n — 1 to size n and the number of 
height-labeled Dyck n-paths is indeed {2n — 1)!!. 

2.7 Height-labeled ordered trees 

Ordered trees of n edges in which each non-root vertex is labeled with a positive integer 
< its height (distance from root). 




1 1 



The "accordion" bijection from ordered trees to Dyck paths — burrow up the branches 
from the root and open out the tree as illustrated — sends non-root vertices to tops of 
upsteps and preserves height. 




The "accordion" bijection 



Thus height-labeled ordered trees transparently correspond to height-labeled Dyck paths. 
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2.8 Overhang paths 

Lattice paths of steps (1, 1), (1, —1), (—1, 1) from (0,0) to (2n, 0) that he in the first quad- 
rant and do not self- intersect [9]. 

Listing the ordinates of the upstep tops is a bijection to trapezoidal words. 

3 Pfaffians, Hafnians and Dyck paths 

The Pfaffian is usually defined for a skew-symmetric matrix but it can just as well be 
defined for the upper triangular array T = (a;ij)i<j<j<2n : 

Pf (T) = ^ Sgn((a, h))Xa{l)b{l)Xa{2)b{2) ■ ■ ■ Xa(n)b{n), 
{a,b)GP„ 

a sum over all (2n — 1)!! perfect matching permutations (a, b) in Vn (where sgn((a, b)) 
is the sign of the permutation). The Hafnian of T is given by the same sum but without 
the signs: 

Hf(T)= 

Xa{l)b{l)Xa{2)b{2) " " ' Xa(n)b{n) 

(a,b)GP„ 

Obviously, the Hafnian of the all Is array is (2n — 1)!!. 

Proposition 1. For an array T = (xjj)i<j<j<2n with constant rows Xtj := Xi, 

n 

where An denotes the set of increasing sequences a = {a{l) , a{2) , . . . , a{n)) satisfying 
1 < a{{} <2i-l. 

Proof. The map (a, b) ^ a maps perfect matching permutations onto An- Thus every 
term in Hf (T) has the form HiLi -^"(O f*-*^ some a G An and the question is, how many 
of each form? Given a G An, set R = {a(l),a(2), . . . ,a{n)} and C = [2n] \ R, and let 
Ci = {c E C : c > a{i)} and Cj = | Cj |. Clearly, C„ C C„_i C . . . C Ci. The fej's in a 
perfect matching permutation (a, b) are subject only to the two restrictions: bi G Ci and 
all 6j's distinct. Thus there are c„ choices for 6„, c„_i — 1 choices for 6„_i, c„_2 — 2 choices 
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for bn-2 and so on. Hence the coefficient of nr=i -^aC*) -^^ (-^) YYi=i {^i — — i)) and 
a simple check shows that Ci — {n — i) = 2i — a{i). □ 

A Dyck path can be coded by the positions in the path of its upsteps and this coding is a 
bijection from Dyck n-paths onto An- A Dyck path can also be coded by the heights of the 
tops of its upsteps, giving a bijection from Dyck ra-paths to the sequences Bn = { 
satisfying 6(1) = 1 and 1 < b{i + 1) < b{i) + 1 for 1 < i < n — 1. The two codes are 
related by the equality b{i) = 2i — a{i) for all i. 

positions of upsteps a = (1, 2, 4, 5, 6, 11) 
heights of upstep tops b = (1, 2, 2, 3, 4, 1) 

a Dyck 6-path 

The correspondences, Dyck n-paths An — > Bn-, along with the assertion of Prop. 1 
and the fact that the Hafnian of the all Is array is {2n — 1)!!, now establish item 2.6. 

This proof can be worked up into a bijection from height-labeled Dyck paths to perfect 
matchings (a, b). The positions of the upsteps in the Dyck path give a, and b is formed 
from the positions of the downsteps, using the sequence of height labels h = {h{i))^__^, as 
follows. First, write the downstep position list d in decreasing order. Then b{n) is the 
h{n)-th entry in d, b{n — 1) is the h{n — l)-th entry in the remaining n — 1 elements of d, 
b{n — 2) is the h{n — 2)-th entry in the remaining n — 2 elements and so on. An example 
is illustrated. 

a= (1,2,5,6,7, 10,12,13), 
d = (16,15,14,11,9,8,4,3), 
h = (1,2, 1,1, 3, 1,2, 2) 

upsteps, downsteps, labels 

So 6(8) is the second entry of d, namely 15; 6(7) is the second entry of (16, 14, 11, . . .), 
namely 14; 6(3) = 16 and so on. The result is b = (4, 3, 9, 11, 8, 16, 14, 15) and the perfect 
matching is 

1 4/2 3/5 9/6 11/7 8/10 16/12 14/13 15. 
There is an analogue of Prop. 1 for the Pfaffian. 





1 2 3 1 5 6 7 s in 11 12 13 14 15 16 

height-labeled Dyck 8-path 
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Proposition 2. For an array T = (a;jj)i<j<_,<2n with constant rows Xij := Xi, 

Pf (T) = X1X3X5 . ..X2n-1- 

Proof. For a perfect matching permutation (a, b), take the smallest i for which a{i) 
and b{i) are not consecutive integers (if there is one). Then the entries in (a, b) up through 
a{i) necessarily form an initial segment of the positive integers ending at a{i) = 2i — 1, 
and a{i + 1) = 2i. So interchanging b{i) and b{i + 1) gives another perfect matching 
permutation. Both contribute the same product to the Pfafiian but with opposite signs 
and hence they cancel out. The only surviving permutation under this involution is the 
identity, which contributes xix^x^ . . . X2n-i- □ 

Corollary 3. For the array T = (xjj)i<i<j<2n with Xij := i, 

Pf(T) = (2n- 1)!!. 

Proof. Put Xi = i in Prop. 2. 



4 Round-Number Identities 
4.1 

T.iklil^^f'-^y-'-i^''-^^- 3)" = {^n-l)U (1) 

This identity counts increasing ordered trees of size n by size k of the leftmost subtree 
of the root. To see this, simply condition on the vertex set of the leftmost subtree. The 
bivariate generating function J2n>i k>o ik+i) (^^ ~ l)!!(2n — 2k — 3)\\j^y'^ is 

1 - y/l - 2xy 
yVl - 2x ' 

and the first few values are 



^\k 
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2 


3 


4 


1 


1 










2 


2 


1 








3 


9 


3 


3 






4 


60 


18 


12 


15 




5 


525 


150 


90 


75 


105 
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4.2 



^ /2n - A; - 1\ (2n - 2k - l)(2n - k + 1] , 

^[ k-i )- rn ^<'" 



2A;-3)!! = {2n-V)\\ 



(2) 



This identity counts increasing ordered trees of size n by length k of the rightmost path 
from the root — the path that starts at the root and successively goes to the rightmost 
child until it reaches a leaf. To see this, let uin^ k) be the number of such trees. Clearly, 
u{0, 0) = 1 and, for n > 1, u{n, 1) = n{2n — 3)U since there are n choices for the rightmost 
child of the root and this child is a leaf. Now suppose n > k > 2 and consider trees of 
size n—1. If the rightmost path has length > k — 1 then there is just one way to insert n 
to get a size-n tree with rightmost path of length k, and n then ends the rightmost path. 
On the other hand, if the rightmost path has length k, then adding n anywhere except 
as the rightmost child of one of the k + 1 vertices on the rightmost path gives a size-n 
tree with rightmost path of length k, and in this case n does not end the rightmost path. 
Thus we have the recurrence 

u{n,l) = n{2n-3)\\ 

n-l 

u{n,k) = u{n — 1, j) + {2n — k — 2)u{n — 1, k) n>k>2, 

j=k-i 

and the summand in (2) satisfies this recurrence. The recurrence leads to a partial differ- 
ential equation for the generating function F{x,y) = X]n>fc>o^('^' ^)"nr^'^' 



(1 -y){l- 2x)F,{x, y) + y(l - y)F,{x, y) + y^F{x, y) = y/VT^, 
with solution 

F{x,y) 

The first few values of u{n, k) are 



1 - (1 -y)ef(i-^T^) 
y^l - 2x 








1 


2 


3 4 5 





1 








1 





1 






2 





2 


1 




3 





9 


5 


1 


4 





60 


35 


9 1 


5 





525 


315 


90 14 1 


k = 


in 


the array 


wants to be 



nice generating function.) 
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Remark The summand in ( ) can be written somewhat more compactly by distin- 
guishing the cases k even or odd: the summand is 

Our interpretation of (j) is equivalent to the assertion that the number of increasing 
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ordered trees of size n whose rightmost path from the root has length > is (^"^ ) (2 
2k — 1)!!. We now prove this assertion bijectively. Form a "vertex" set V = {1, 2, . . . ,n} 
and an "edge" set E = {1,2, . . . ,n — k} and make them disjoint by using subscripts 
V and E on their entries. Thus \ V U E \ = 2n — k. It suffices to exhibit a bijection 
from the trees being counted to pairs consisting of a fc-element subset X oi V U E and a 
standard increasing ordered tree Tq of size n — k, since these pairs are clearly counted by 
(^"'~^)(2n — 2k — 1)!!. We use the tree illustrated in Fig. la) as a working example with 
n = 14 and k = 5. 



12 9 10 12 9 10 






Fig. la) Fig. lb) 



Call the first k edges of the rightmost path from the root the base edges and their child 
vertices the base vertices (2,4,6,7, 10). Say a base vertex is fertile if it has a child that 
is not a base vertex and barren otherwise. Highlight each barren vertex and the leftmost 
edge from each fertile vertex, and color the base edges red as in Fig lb). Thus there are 
k highlighted vertices/edges. The labels on the highlighted vertices are themselves the 
contribution from V to the fc-element set X. The red edges will be deleted by an iterative 
cut-and-paste procedure to get the required (n — k)-edge tree Tq and then the positions 
of the highlighted edges in Tq will determine the contribution from E to X. 

First, erase the barren vertices and standardize the vertex labeling — replace smallest 
by 0, next smallest by 1, and so on, to get the first tree in Fig. 2. 
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Fig. 2 

Now eliminate the remaining red edges, one at a time, in increasing order of their 
child vertices. Let b denote the current smallest base vertex (initially b = 2). Cut out the 
subtree rooted at b, discarding the label b and its parent edge, and re-root this subtree 
at b — 1, placed so that it lies to the right of the existing subtree rooted at b — 1, and 
standardize vertex labels. Repeat this process on the second, third, . . . , red edge in 
turn to obtain the desired (n — fc)-edge tree Tq (the last tree in Fig. 2) with some edges 
highlighted. Observe that the progress of individual highlighted edges in the successive 
trees of Fig. 2 can be discerned even though the labels on their endpoints will change. 

Finally, list the edges of Tq in standard order, that is, in increasing order of their 
parent vertices (preserving, of course, the order of edges with a common parent vertex), 
and take the positions of the highlighted edges in this list as the contribution from E to 
the /c-element set X. In the example, Tq has defining edge list 

((0, 3), (0, 1), (1, 9), (1, 2), (3, 8), (3, 4), (3, 6), (3, 5), (4, 7)) 

and the highlighted edges — (1, 9), (3, 4), (3, 5) — are in positions 3,6,8. The net result is 
the pair {X, Tq) with X = {4y, lOy, 3^;, 6^;, 8^;} and Tq as just given. 

Is this process reversible? The entries in X with subscript E determine the highlighted 
edges in Tq, and the entire process can now be reversed step- by-step provided we know the 
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order in which the highhghted edges were originally processed. But this order is precisely 
their left-to- right order in the standard hsting of the edges of To- Q 

Remark (^"'^'^)(2n — 2k — 1)!! is also the number of increasing ordered trees of size 
n + 1 whose root has k + 1 children, increasing from left to right: given an increasing 
ordered tree of size n whose rightmost path from the root has length > k, delete the base 
edges, increment all labels by 1, and then attach the old root and the k base vertices in 
increasing order to a new root 0. This is a bijection to the trees in question. 



4.3 



n / _ 1 \1 

Ef^^(2^-3)!! = (2n-l)!! 



(3) 



This identity counts increasing ordered trees of size n by the maximum k of the young 
leaves where a young leaf is a leaf with no left sibling (leaf or otherwise). Every nonempty 
increasing ordered tree has at least one young leaf. To establish this count, consider the 
2n — 1 ways to produce an increasing ordered tree of n edges by inserting a leaf n into an 
increasing ordered tree of n — 1 edges: either as the leftmost child of one of the n vertices 
or so that the new edge is the immediate right neighbor edge of one of the n — 1 existing 
edges. In the first case, the maximum young leaf becomes n; in the second, it is preserved. 
In particular, the maximum young leaf becomes n in n ways from any increasing ordered 
tree of n — 1 edges. Hence, if u{n, k) denotes the number of increasing ordered trees 
with n edges whose maximum young leaf is A;, we have u{n,n) = n{2n — 3)!! and the 
recurrence u{n, k) = (n — l)u{n — 1, /c) for 1 < /c < n. Iterating the recurrence yields 
u{n, k) = {n - l)'-^k{2k - 3)!! for all 1 < A; < n. 

The bivariate generating function Ylm k>i _ jji ^ (2A; — 3)!! 



and the first few values are 
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x){l 


-2X2/) 3/2' 
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3 4 5 
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1 
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2 
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12 


27 60 
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24 


48 


108 240 525 
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4.4 



n/2 
k=0 



n 
2k 



2k 
k 



^ = (2^-1)" 



(4) 



This identity counts perfect matchings of \2n] by number k of matches in which both 
entries are < n. In fact, (4) is the special case r = of a family of identities indexed by 
a nonnegative integer r: 



(n+r)/2 

E 

k=r 



n 

2k -r 



2k -r 
k 



{n + r)- (n — r)!- 



1 



l2fc-r 



(2n- 1)!! 



(5) 



where [n + r)- = + r)(n + r — 1' 
a straightforward direct count shows 



{n + 1) is the falling factorial. For given r, as 
) counts perfect matchings of [2n] by number k 
of matches in which both entries are < n + r (also by number k of matches in which 
both entries are < n — r). The bivariate generating function Ylm k>r (^VO ~ 



(2r- 1)!! 



((1 



(2r+l)/2 • 



For r = 0, the first few values are 



^\k 
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3 





1 
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1 
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1 






3 


6 


9 






4 


24 


72 


9 




5 


120 


600 


225 




6 


720 


5400 


4050 


225 



and for r = 1, the first few values are 
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2 


3 


1 


1 






2 


3 






3 


12 


3 




4 


60 


45 




5 


360 


540 


45 


6 


2520 


6300 


1575 
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By transposing factors, ( ) can be written in the alternative form 




2k — r ) \ k 




(6) 



an identity that counts bicolored UDF paths of length n ending at height r by number 
k of upsteps. A UDF path is a lattice path of upsteps (1, 1), downsteps (1,-1), and 
flatsteps (1,0); bicolored means each fiatstep is colored red or blue; height is relative to 
the horizontal line through the initial vertex. 

The case r = 1 of ( ) also counts increasing ordered trees of n edges by number k of 
young leaves (as defined in Section 1.3). To see this, note that inserting a leaf n into an 
increasing ordered tree of n — 1 edges {2n — 1 possible ways) always either preserves or 
increments (by 1) the number of young leaves, and the number of ways to increment is 
n — 2k where k is the number of young leaves. This observation leads to the recurrence 

u{n, k) = {n + 2k- l)u{n -l,k) + {n-2k + 2)u{n - 1, - 1) I <k <{n + l)/2 

for the number u{n, k) of increasing ordered trees of n edges and k young leaves, and the 
recurrence is satisfied by the summand. 



This identity counts 

(1) Stirling permutations of size n by first entry k, 

(2) Stirling permutations of size n by position 2A; — 1 of the first 1 (the position of the 
first 1 is necessarily odd), 

(3) increasing ordered trees of size n by the parent — 1 of n, 

(4) increasing ordered trees of size n by the label k on the leaf in the minimal path from 
the root. The minimal path starts at the root and successively travels to the smallest 
child vertex until it arrives at a leaf. 



(1) and (4) The number of choices in building up the object, inserting a pair ii in the 
permutation or a leaf i in the tree for i = 1 to n, is successively 1, 3, . . . , 2/c — 3, 1, 2k, 2k + 
2, . . . , 2ra — 2 and their product is the summand. 



4.5 




(7) 



k=l 



Proofs 
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(2) If a letter occurs to the left of the first 1, then both occurrences do so, and so the 
number of such permutations is (^ij) [choose support set for the first 2k — 2 entries] 
x(2/c — 3)!! [form a Stirling permutation on this support] x(2n — 2k)\\ [form a Stirling 
permutation of size n — k + 1 that starts with a 1], and {^^^{2k — 3)!!(2ri — 2k)\\ is an 
equivalent expression for the summand. 

(3) Here, build up the tree by successively inserting leaves 1, 2, . . . , /c — 1, then n must 
be inserted as a child of A; — 1, then proceed to insert /c. A; + 1, . . . , n — 1. The number of 
choices is the same as in (1) and (4) above. 

The bivariate generating function Y.n>i,k>i ^^"'(2fc-2)!7 (rf- 1)! ^''"^ 

1 

(1 -2x)^r^^2i^' 

and the first few values are 



^\k 


1 


2 


3 


4 


5 


1 


1 










2 


2 


1 








3 


8 


4 


3 






4 


48 


24 


18 


15 




5 


384 


192 


144 


120 


105 



With rows reversed, this array is entry A122774 in OEIS [l(J], and the reversed array 
counts (i) height-labeled Dyck paths by the position among the upsteps of the last upstep 
with label 1, and (ii) increasing ordered trees by the maximum child of 1. 



This identity counts increasing ordered trees of size n by smallest child k of 1 {k = 
1 if vertex 1 has no children). To see this, let u{n,k) be the number of such trees. 
Consideration of the effect of inserting a leaf n into an increasing ordered tree of size n — 1 
on the smallest child of 1 yields the recurrence 

u{n,l) = {2n-2)\\ 

u{n,k) = {2n-l)u{n-l,k) 2<k<n-l 
u{n,n) = (2n-4)!!, 
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satisfied by the summands. The recurrence leads to a first-order ordinary differential 

equation for the generating function F{x,y) — X^„>fc>iM('^, k) , \ . y'^~^: 

- - [n i). 



(1 - 2x)F^ - 3F + 1/(1 - 2x) 



1 — 2xy ' 



with solution 

F{x,y) = 



(l-2x)3/2 (l-2x)3/2y 

The first few values of u{n, k) are 



^ ^y-y^ ^^^-1 / ^y-y'' (-i + 2xy + ^/i^) 



l-2y + 2xy^ 
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2 


3 


4 


5 


1 


1 










2 


2 


1 








3 


8 


5 


2 






4 


48 


35 


14 


8 




5 


384 


315 


126 


72 


48 



With rows reversed, this array has an interpretation for which the recurrence relation 
is not immediately obvious: let v{n, k) denote the number of Stirling permutations of size 
n for which the maximum M of the entries preceding the first 1 (taken as n + 1 if the 
permutation starts with a 1) is A; (2 < A; < n + 1). Then v{n, k) — u{n, n + 2 — k). This 
follows from the next Proposition by comparing recurrences. 



Proposition 4. 



(i) v{n,n + l) 
(ii) v{n, 2) 
{in) v{n, k) 



(2n-2)!!, 
(2n-4)!!, and 

{2n-l)v{n-l,k-l) for3<k<n. 



Proof 

(i) Consider permutations that start with a 1. Inserting two adjacent n's immediately 
after any one of the 2n — 2 entries in such a permutation of size n — 1 gives one of size n. 
Thus v{n, n + 1) — {2n — 2)v{n — l,n), and (i) follows. 

(ii) A Stirhng permutation with M — 2 necessarily begins 221 Deleting the initial 

2s is a bijection to Stirling permutations of size n — 1 that start with a 1, counted by 
v{n-l,n) = (2n-4)!!. 
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(iii) Now suppose 3 < A; < n. A Stirling permutation of size n—1 with M — k — 1 yields 
2n — 1 Stirling permutations of size n with M — k (all distinct) as follows. Increment 
each entry > 2 by 1 and tentatively insert a pair of adjacent 2s anywhere in the resulting 
permutation (2n — 1 choices) to obtain a permutation a. Let / denote the initial segment 
of (7 terminating at the first 2. If no entry exceeding 2 occurs exactly once in I, then a 
is already a Stirling permutation and we may leave the tentative 2s in place. Otherwise, 
choose the smallest i > 2 that occurs exactly once in /, say ii, and interchange the 2s and 
iis in a to obtain a perm ui. Let Ii denote the initial segment of ai terminating at the first 
ii. If no entry exceeding ii occurs exactly once in Ii, then ai is a Stirling permutation, 
and stop. Otherwise, proceed similarly to get i2 > ii, interchange the iis and to get 
(72 and continue until you arrive at a permutation ak that does have the Stirling property. 
The original Stirling permutation of size n — 1 and the location of the tentative 2s can be 
recovered from ak and so this process is a bijection V{n — 1, k — 1) x [2n — 1] ^ V(n, k) 
where V{n, k) is the set counted by k). □ 

Also, v{n, k) is the number of increasing ordered trees of size n with k the maximal 
descendant of 1 (taken as n+1 if 1 is a leaf since "descendant" here means "proper descen- 
dant"). The generating function Y^n>i X]fc=2 '"(^^ ^) (^^1 1)1 v{n,k) is marginally 
more concise than that for u{n, k): 

(1 - 2a;|/)3/2 (1 - 2xy)y^ ^ \ 2 - 2x - y )' 



4.7 

p»-3)"+E (,,|i)(J:i) -p»-i)" (9) 

This identity counts Stirling permutations of size n by smallest entry k following the last 
n {k — Q ii the last entry is n) . The recurrence relation for these permutations is 

ii(n,0) = (2n-3)!! 

u{n,k) = {2n - l)u{n - l,k) l<k<n-2 
u{n,n-l) = 2(2n-5)!!, 

satisfied by the summand. 

Remark The identity itself is trivial to prove, since the sum is telescoping, but the 
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generating function J2n>i 



, 0<k<n 



^-^ u{n, k) 



X" 



(n-1) 



■y'' is cumbersome: 



1 



+ 



I yi - 2xy 

;i-2x)3/2 ' VI -2x l-2x ' 

1-2/ 



2^(1 - 2s)3/2 
The first few values of u{n, k) are 



log 



l + y-4xy- 2^(1 - 2x)y{l - 2xy) 

(1 - Vy)' 



^\k 





1 


2 


3 


4 


1 


1 










2 


1 


2 








3 


3 


10 


2 






4 


15 


70 


14 


6 




5 


105 


630 


126 


54 


30 



4.8 



fc=i 



2n-k 

k-: 



{2n -2k- 1)!! = (2n - 1) 



(10) 



This identity counts increasing ordered trees of n edges by outdegree k of the root. 
The k\ factor allows us to assume that the k children of the root increase from left to 
right, and the number of such trees is given by the Remark at the end of Section 4.2. 

The identity also counts height-labeled Dyck n-paths by length k of first ascent. To 
show this, we exhibit a bijection from height-labeled Dyck n-paths to increasing ordered 
77,-trees that sends "length of first ascent" to "number of children of the root" . The tree 
is built up using trees with some of their leaves unlabeled. The construction involves a 
pair of sequences (a(z))"__^ and that characterizes the path: 

a{i) = #upsteps immediately preceding the ith downstep, and 
b{i) = label on the upstep matching the ith downstep. 

Thus < a{i) < n and Yli=i '^(^) — ^- -^^ ^^ep i {1 < i < n), attach a{i) unlabeled leaves 
to (leaf) vertex i — 1 and apply label i to the b{i)-th unlabeled leaf in walkaround order 
as illustrated. 




i 


1 


2 


3 


4 


5 


6 


7 


a{i) 


2 


3 





1 








1 




2 


3 


3 


2 


1 


1 


1 



height-labeled Dyck 7-path 
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defining sequences 






step 1 









construction of corresponding tree 



We remark without proof that the identity also counts leaf-labeled 0-2 trees of size n 
by length of the path from 1 to the root. 



The bivariate generating function En>o,ik>o^'fT-i ^) (2^ - 2A; - l)n|Jy'= i 

1 - 1/ - y yi - 2a; 
l-2y + 2xy^ 



IS 



Omitting the empty tree {n = k = 0), the generating function Yln>i k>i "fe-i ) (^^ ~ 

2k — , im V^"^ looks rather different: 
' (n— 1)!^ 



and the first few values are 



VI- 


2x (1 


-y 


^y^\-2xy' 




1 


2 


3 4 5 


1 


1 






2 


1 


2 




3 


3 


6 


6 


4 


15 


30 


36 24 


5 


105 


210 


270 240 120 
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This array is entry A102625 in OEIS. 



The result of Section I S can be refined to give the joint distribution of first ascent 
length and first peak upstep label: let u{n,j,k) denote the number of height-labeled 
Dyck ra-paths whose first peak upstep has label j and whose initial ascent has length 
k {n > k > j > 1). Then, since j is uniformly distributed over [k], u{n,j,k) = 
[k — {2n — 2k — 1)!! for 1 < j < /c, leading to the generating function 

1 

VI - 2x {l-yz + yzVl - 2x) [l - z + z^\ - 2x) ' 

The bijection of this section shows that u{n,j,k) also counts increasing ordered trees 
whose root has k children among which vertex 1 is the j-th. 



5 Non-Round Identities 
5.1 

n\ + J2 ((^ - 1)"(2^ - - k\\{2n -k- 1)!!) = {2n - 1)!! (11) 

k=l 

This identity is trivial to prove — the sum is telescoping and collapses to {2n — 1)!! — 
{n — l)!!n!! — but it has an interesting interpretation: it counts height-labeled Dyck n- 
paths by length k of the first descent where a descent is a maximal sequence of contiguous 
downsteps and the term n! corresponds to k = n. To see this, let u{n, k) denote the number 
of height-labeled Dyck n-paths whose first descent has length >k{l<k<n). The next- 
size-up construction described in fi) yields one of these paths precisely when the (n— 1)- 
path has first descent of length > k and the specified vertex is not the terminal vertex of 
one of the first k — 1 downsteps in the first descent. Thus u(n, k) = (^{2n—l)~{k — l))u{n — 
1, k) = {2n — k)u{n — k) for n > k. Together with the obvious initial case u{k, k) = k\, 
this recurrence yields that u{n, k) = k\{k + 2)""'^'^, where fc"'^ = k{k + 2){k + A) ■ ■ ■ to n 
factors is the rising double factorial. Equivalently, u{n, k) = {k — l)!!(2n — k)\\. Thus the 
number of height-labeled Dyck n-paths with first descent of length k is u{n,n) = n\ for 
k = n, andM(n, k)-u{n,k + l) = (A; - l)!!(2n - A;)!! - A;!!(2n - A; - 1)!! for 1 < A; < n-1. 
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and the first few values are 



r) \ ^ 


1 
1 


9 


Q 
O 


A 


o 


1 


1 










2 


1 


2 








3 


7 


2 


6 






4 


57 


18 


6 


24 




5 


561 


174 


66 


24 


120 



Here <[^]> is tlie second-order Eulerian number (indexed so tliat 1 < A; < n) A008517. 
Tfiis identity counts 

(1) Stirling permutations of size n by number k of descents (including a conventional 
descent at the end), 

(2) Stirling permutations of size n by number k of plateaus, that is, pairs of adjacent 
equal entries, 

(3) increasing ordered trees of size n by number k of leaves, 

(4) height-labeled Dyck n-paths by number k of upstep-free vertices, and 

(5) trapezoidal words of length n by number k of distinct entries. 
Furthermore, the second-order Eulerian triangle with reversed rows counts 

(6) increasing ordered trees by number of descents where a descent in a tree is a pair of 
adjacent sibling vertices with the first larger than the second (no conventional descents), 
and 

(7) height-labeled Dyck ra-paths by number of peaks. 
Proofs 

(1) Several proofs are given in [2]. 

(2) Stirling permutations are usually defined on support set [n] but of course can be sim- 
ilarly defined on an arbitrary set of positive integers. We present a bijection 0, actually 
an involution, on Stirling permutations of arbitrary support set that preserves size and 
interchanges descents" and plateaus". First, (p is the identity on the empty permu- 
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tation. A nonempty Stirling permutation can be written as the concatenation AmBmC 
where m is the smallest integer in its support set and A, B, C are perforce themselves 
Stirling permutations. Now define recursively by 

(f){AmBmC) = 0(A) m 0(C) m 0(5). 

One checks, using induction, that has the required properties. 

(3) This result follows from Janson's bijection (see Section ') since it clearly sends leaves 
to plateaus. 

(4) There are 2n+l vertices in a Dyck n-path and a vertex is upstep-free if it is not incident 
with an upstep. The recurrence of {-.{)) for height-labeled Dyck paths can easily be refined 
to a recurrence for h{n, k), the number of height-labeled Dyck n-paths with k upstep-free 
vertices, which turns out to be h{n,k) = kh{n — l,k) + {2n — k)h{n — l,k — 1) — the 
defining recurrence for second-order Eulerian numbers [2, p. 27]. 

The recurrence leads to a bijection from height-labeled Dyck n-paths to Stirling per- 
mutations that sends # upstep-free vertices to # conventional descents. The bijection uses 
identical "codings" for the two classes of objects. Consider the alphabet consisting of two 
copies of the positive integers, distinguished by subscripts Y and A^. Let denote the 
set of words (w{i))^_^ over this alphabet satisfying the following condition for each i. If 
w{i) has subscript Y, then its value is < 1 + # subscripts preceding w{i), while if w{i) 
has subscript A^, then its value is < 2i — 2 — # subscripts A^ preceding w{i). Clearly, there 
are 2i — 1 choices for w{i) regardless of the preceding entries, and so | C„ | = {2n — 1)!!. 
Indeed, a word w in Cn can be represented as a lattice path along with a set of lattice 
points: discard u'(l) — necessarily ly — and subtract 1 from each remaining value to get 
a sequence of n — 1 nonnegative integers [b{i)) and a sequence of n — 1 subscripts. The 
subscripts give the path letting Y denote a downstep (1, —1) and A^ a fiatstep (1, 0), and 
the ith lattice point is b{i) units above (resp. below) the terminal point of the ith step if 
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the corresponding subscript is (resp. Y), as illustrated. 



(ly Iat 2y Sat 2^ 3y) 



i 




N Y N N Y 



14 12 




word in C, 



•n 



path and set of lattice points 



The ordinates of the lattice points with a prepended form a symmetric trapezoidal 
word (Section 2.1). The path is redundant since it can be recovered from the lattice points 
and so we have an explicit bijection from C„ to symmetric trapezoidal words. 

To code a height-labeled Dyck ra-path by a word in C„, delete the last upstep U with 
label 1 along with the last downstep and decrement labels after U to get an {n — 1)- 
path Pi with a distinguished vertex v (where the upstep was deleted). If v is upstep-free 
(resp. upstep-incident) in Pi, the subscript on w{n) is Y (resp. A^) and its value is the 
number of upstep-free (resp. upstep-incident) vertices weakly preceding v. Repeat on Pi 
to get w{n — 1) and so on, ending with w{l) := ly. This gives a bijection that sends 
# upstep-free vertices to # subscripts Y. 

Likewise, to code a Stirling ra-permutation by a word in C„, delete the two ns (nec- 
essarily adjacent) to get an {n — l)-permutation with a distinguished gap (possibly at 
either end). Recalling that the last gap is considered a descent, if the distinguished gap 
is a descent (resp. non-descent), the subscript on w{n) is Y (resp. A^) and its value is 
the number of descent (resp. non-descent) gaps weakly preceding the distinguished gap. 
Repeat to get w{n — 1) and so on, again ending with w{l) := ly. This gives a bijection 
that sends # descents to # subscripts Y. 

(5) This result is stated without proof in [ ]. A Stirling permutation can be built up in a 
unique way by starting with a plateau of two Is, inserting a plateau 22 in one of the 3 gaps 
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in -1-1-, then inserting a plateau 33 in one of 5 gaps and so on. Define a mapping from 
Stirling permutations of size n to trapezoidal words w — (tyi)F=i ^ follows. Set wi — 1. 
If 22 is placed between the two Is, set ^2 = 1; if 22 is placed to the left of the ones, set 
W2 — 2 (the smallest number not yet appearing in w), else W2 — S. In general, if kk is 
placed inside a plateau, say the ith plateau (left to right), Wk is the ith smallest number 
already appearing in w; otherwise kk is placed in one of the remaining gaps, say in the 
jth of the remaining gaps, and Wk is the jth smallest positive integer not yet appearing 
in w. For example, 5512234431 11442 as follows. 

11 1221 122331 ^ 12234431 ^ 5512234431 

Wi — 1 W2 — 1 Ws — A* W4 — — 2 

* 4 is the third smallest number not yet appearing in w 
^ 4 is the second smallest number already appearing in w 

It is easy to check that this algorithm defines a bijection from Stirling permutations of 
size n to trapezoidal words of length n and that it sends plateaus" in the permutation 
to distinct entries" in the word. 

(6) Translated to Stirling permutations using Janson's bijection described above, descents 
in a tree become strong descents in the permutation where a strong descent is a descent 
a > b involving the first of the two occurrences of b in the permutation. Note that the 
analogous notion of strong ascent is superfluous because all ascents in a Stirling permu- 
tation are strong. Consideration of the effect of inserting nn into a Stirling permutation 
of size n — 1 on the number of strong descents leads to the defining recurrence for the 
reversed second-order Eulerian triangle using the following fact, proved by induction: in 
a Stirling permutation of size n, # strong descents -|- # ascents —n — 1. 

(7) The number of upstep-free vertices in a Dyck n-path is related to the number of peaks: 
their sum is 1. It follows that the number of height-labeled Dyck n-paths with k peaks 
is the second-order Eulerian triangle with reversed rows. 



2"-*^ = (2n-l)!! (13) 

Here [^] is the Stirling cycle number (and [q] = 1). This identity counts 

(1) Stirling permutations of size n by number k of left-to- right (LR) minima. (For example, 
5523443211 has 3 LR minima, namely 5,2,1), 
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(2) increasing ordered trees of size n by number k of edges in the minimal path from the 
root as defined in Section 4.5, and 

(3) height-labeled Dyck n-paths by number k of upsteps in the first ascent with label 1. 
Proofs 

(1) It suffices to exhibit a bijection from Stirling permutations a of size n with k left- 
to-right minima to pairs {A, r) where A is a subset of [n — k] and r is a permutation 
of [n] with k left-to-right minima (recall the number of such r is [^]). Split a just 
before each LR minimum to write a as a concatenation ai, cr2, . . . , 0"^ of lists. Each ai 
is itself a Stirling permutation and their support sets, taken in order, form a partition 
of [n] into blocks whose smallest elements are decreasing left to right. In view of these 
observations, the general case follows from the special case = 1 by amalgamating subsets 
and concatenating permutations, both taken on the appropriate support sets. So suppose 
/c = 1, implying that a starts with a 1. For j G [2, n], let ij denote the last entry preceding 
j in a that is < j. Thus, for a = 12 5 5 2144 3 3, 12,^3,^4 all = 1 and is = 2. Now take 
A = {j & [2, n] : both occurrences of ij in a precede the first occurrence of j}, a subset of 
the (n — l)-element support set [2,n]. As for r, observe that {ij)^=2 ^ [1] x [2] x ■ ■ • x [n — 1] 
and so corresponds to a permutation of the (n — l)-element set [2, n]. Prepend a 1 to this 
permutation to get the permutation r of [n] with k = 1 LR minima. We leave the reader 
to verify that a can be recovered from the pair {A, r). 

(2) Adding a leaf n to a tree of size n — 1 preserves the length of the minimal path except 
when n is added as a child of the leaf that terminates the minimal path. We thus have 
the recurrence 

u{n, k) = {2n — 2)u{n — 1, A;) + u{n — l,k — 1), 

satisfied by the summand because it reduces to the basic recurrence for the Stirling cycle 
numbers: K] = (n - 1) M + K] • 

(3) We give yet another bijection from height-labeled Dyck paths to increasing ordered 
trees. This one sends # upsteps immediately preceding the ith downstep to the outdegree 
of vertex i — 1 (1 < i < n) — in brief, full ascent sequence of path —>■ fertility sequence of 
tree — and sends locations of Is on first ascent to locations of LR minima among children 
of the root. In particular, 7^ Is = #LR minima, and Janson's bijection identifies LR 
minima among children of the root with left-to-right minima in a Stirling permutation. 
So (3) will follow from (1). The bijection is an algorithm to generate the edge list. With 
n = 10, and the path below working example. 
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height-labeled Dyck 10-path 



list the total number of downsteps preceding the i-th upstep for z = n, n — 1, . . . , 1 (column 
2 in Fig. 3) and the label on the 2-th upstep (column 3). 



i 


#Ds 
< i-th U 


label on 
i-th U 


candidate 
children 


selected 
edge 


10 


7 


3 


8 910 


710 


9 


7 


1 


89 


78 


8 


3 


5 


45679 


39 


7 


3 


1 


4567 


34 


6 


3 


2 


567 


36 


5 


2 


3 


357 


27 


4 





2 


1235 


02 


3 





1 


135 


01 


2 





2 


35 


05 


1 





1 


3 


03 




Fig. 3 



Each entry j in column 2 is the parent vertex of an edge. The candidates for its child 
are the entries of [j + l,n] not already having a parent. Initially j = 8, no vertex has a 
parent, and the candidate children are 8,9,10. The corresponding label in column 3 then 
determines the child, here label 3 selects the third candidate, namely 10, so 7 10 becomes 
an edge; 10 now has a parent and is deleted from the candidate set. The label 1 in the next 
row selects the first candidate, namely 8, and adds 7 8 to the edge list. Proceed similarly 
to get all n edges. Then the last column, read upwards, is the edge list in standard order, 
giving the tree shown. 

The generating function ^22=1 [^]2"'~^x"' for row n is nr=o^('^ + bivariate 
generating function E„,yt>o 2""^lr2/'' (1 - 2xyy/^. See A039683 in OEIS. 
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Consider the following modification of the minimal path from the root in increasing or- 
dered trees: the right-then-minimal path is the path that goes from the root to its right- 
most child, then follows minimal children to a leaf. Let u{n, k) denote the number of 
increasing ordered trees of size n whose right-then- minimal path has length k, so that 



^u{n,k) = {2n - 1)!!. 



k=l 



Then u{n, k) satisfies the recurrence 
u{n,l) = n(2n-3)!! 

u{n, k) = {2n — 3)u{n — 1, fc) + u{n — l,k — 1) 2 < k < n, 



leading to a differential equation for the generating function y) = X]n>fc>i '^(''^^ 



2x)F,{x,y) 



with solution 



F{x,y) 

The first few values oi u{n,k) are 



VI -2a; 

y {1 - {1 - 2xy-y/^) 
{2-y)VT^^W 



{l-y)F{x,y), 



^\k 


1 


2 


3 


4 5 


1 


1 








2 


2 


1 






3 


9 


5 


1 




4 


60 


34 


10 


1 


5 


525 


298 


104 


17 1 



6 Refinements 
6.1 

The interpretations of Sections 1.3 and 5.3 have a common refinement 

(^n-l)^= (2n- 1)!!. 

The summand is the number of increasing ordered trees of size n whose minimal path 
from the root ends at the leaf labeled k and whose length is j: the number u{n, k,j) of 



k-1 
.J - 1. 
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such paths satisfies the recurrence 



u{n,n,l) = {2n-2)\\, 

u{n,k,l) = k>2, 

u{n,k,j) = {2n — 2)u{n — 1, k, j) 2<j<k<n — 1, 



n-l 



u{n-l,i,j-l) 2<j< 



n. 



i=j-i 



The generating function Ei<j<fc<n"(^' ^'-^^ (n - 1)! ^^ ^ 

1 



F{x,y,z) 



(1 -2x)(l -2xyy/^' 



Likewise, the interpretations of Sections 4.8 and 5.3 have a common refinement: 



l<j<fc<n 



2n-k-l 
k-1 



{2n-2k- 1)!! = (2n - 1)!!. 



The summand f (n, fc, j) is the number of height-labeled Dyck ra-paths whose first ascent 
has length k and contains j Is. The generating function X]o<j<fc<n ^(^' ^i 3)^y^ 

1 - y - yVl - 2x 
1 -2y + 2xy^ 



6.2 

The result of Section 1.8 can be refined to give the joint distribution of first ascent length 
and first peak upstep label: let u{n, j, k) denote the number of height-labeled Dyck n-paths 
whose first peak upstep has label j and whose initial ascent has length k {n > k > j > 1). 
Then, since j is uniformly distributed over [A;], u{n,j, k) = {k — {2n — 2k — l)\\ 

for 1 < j < k, leading to the generating function J2n>ij>i k>i ^(^'i^ ^-^ (rf^T)!^^~^^'^~^ ~ 

1 

VI - 2x {l-yz + yzVl - 2x) {l - z + z^\ - 2x) ' 



6.3 

To refine the result of Section 5.1 and deduce further identities, let u{n,k,j) denote the 
number of height-labeled Dyck n-paths whose first ascent has length = j and first descent 
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has length > k. Thus u{n^ k) = J2]=k ""(^j ^jJ)- Since the first peak in such a path has j 
possible labels, deleting this peak and its label shows that 

u{n, k,j) = ju{n — 1, k — 1, j — 1) for 2 < k < j < n (14) 

In particular, u{n,2,j) is j times the number of size-(n — 1) height-labeled Dyck paths 
with first ascent of length j — 1 and no restriction on the first descent. Hence, by the 
second interpretation of (10), u{'n, 2, j) = j(j — l)(2n — 2 — j)^^^{2n — 2j — 1)!!. This base 
case, together with ( ), yields 

u{n, k,j) = f{2n -k- jy^{2n - 2j - 1)!!. 

Equating '^^=kU{n, k,j) and u{n, k) yields the identity 

n 

^]^{2n-k- j)i^{2n-2j - 1)!! = (A; - l)!!(2n - A;)!!. 

j=k 

Two alternative forms of this identity, eliminating the double factorials, can be found by 
considering the cases where k is even or odd separately: 

with k := 2m, n replaced by n + 2m, and 



^ ,j + 2m + l\^^^j2n-j\ _ m!(2n + 2m + 2): 
-^V j ) n!(2m + l)!(n + m + 1)! 



with k := 2m + 1, replaced by n + 2m + 1. 

This last identity is interesting because it provides another solution to Ira Gessel's 
1987 Monthly Problem E3107: 

Show that is an integer for nonnegative integers m, n. 

Replace m by m — 1 in ( ' ' -), rearrange terms, and cancel a 2 to get 

f] + 2m - 1\ 2j f2n - j\ _ m!(2n + 2m)! 
j ) \fT'~jJ n!(2m)!(n + m)! ' 

exhibiting Gessel's expression as a sum of integers, and this sum is different from that in 
A. A. Jagers' solution [11]. The case m = of ( ) has a simple combinatorial interpreta- 
tion: it counts lattice paths of upsteps and downsteps of length 2n by number of "returns 
to ground level" . 
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